Generalized gauge field theories with non-topological soliton solutions 
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We perform a systematic analysis of the conditions under which generalized gauge field theories of 
compact semisimple Lie groups exhibit electrostatic spherically symmetric non-topological soliton 
solutions in three space dimensions. By the term generalized, we mean that the dynamics of the 
concerned fields is governed by lagrangian densities which are general functions of the quadratic 
field invariants, leading to physically consistent models. The analysis defines exhaustively the class 
of this kind of lagrangian models supporting those soliton solutions and leads to methods for their 
explicit determination. The necessary and sufficient conditions for the linear stability of the finite- 
energy solutions against charge-preserving perturbations are established, going beyond the usual 
Derrick-like criteria, which only provides necessary conditions. 
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The non-linear Born-Infeld (BI) electrodynamics [l| is 
one of the first consistent relativistic field models found 
exhibiting (non-topological) soliton solutions in three 
space dimensions (finite-energy and stable electrostatic 
central-field solutions, as well as dyon-like solutions and 
bi-dyon field configurations 2]). In the last two decades 
there has been a revival of interest on this model and 
its extensions to other fields. These extensions amount 
to obtain new non-linear lagrangian models for a given 
field, by generalizing the (linear or not) original ones 
through the same procedure relating the Maxwell and 
Born-Infeld lagrangians for electrodynamics (involving 
always the well known square root dependence). The 
models so obtained have found applications in many dif- 
ferent physical contexts. We mention here, as examples, 
the generalization of BI procedure to non-abelian gauge 
fields in many space dimensions, suggested by string the- 
ory Q , the analysis of glueballs in terms of soliton so- 
lutions of non-abelian BI gauge field models 0], the use 
of these gauge models for the description of dark energy 
in the context of Cosmology [5], the phenomenological 
analysis of the nucleon structure in terms of BI general- 
ization of the quadratic Skyrme action Q and the search 
for self-gravitating gauge field solutions in General Rela- 
tivity [7|. 

However, the choice of Born and Infold in generaliz- 
ing the action of electrodynamics (as well as the above 
mentioned extensions) which leads to models supporting 
stable finite-energy solutions, is by no means unique. In 
fact the aim of the present work is to determine the class 
of general gauge-invariant lagrangian models (i.e. defined 
as arbitrary functions of the field invariants) which are 
physically consistent and exhibit this kind of soliton so- 
lutions. This is motivated by their potential utility in 
the aforementioned contexts and many others. In Q, 
we have approached this question for the restricted 
cases of (one and many-components) scalar fields a (a; M ), 
with lagrangian densities defined as functions f(X) of 



the kinetic terms A = ^ a d M </> a .<9 M (/> Q . At first glance 
this choice seems rather arbitrary from a physical point 
of view. Nevertheless it is the natural restriction to the 
scalar fields of the general problem considered here for 
the (abelian and non-abelian) gauge fields and, as we 
shall see, those results are necessary for the present anal- 
ysis. Let us summarize here the more useful of those 
results. We have found in Ref.[|| the general conditions 
to be satisfied by "admissible" lagrangian densities f(X) 
determining exhaustively the class of models support- 
ing finite-energy, static spherically symmetric (SSS) so- 
lutions which are linearly stable. We called "admissible" 
the lagrangian densities defined and regular everywhere 
(class- 1 models) or in an open and connected domain in- 
cluding the origin (class-2 models) and leading in both 
cases to energy densities which are positive in these re- 
gions and vanish in vacuum. The SSS solutions obey to 
first integrals of the field equations which only contain 
the field strengths and read 



V„(r)/ 



(1) 



where </> a (r) = f(X) = j| and the A' a s are the 

integration constants, identified as the scalar charges of 
the solutions. The SSS solutions of the many-components 
scalar fields may be explicitly written in terms of those of 
the one-component scalar models, corresponding to the 
same functional form of the lagrangian density, as 



<t>a{r) = ^V(r,A) + Xa, 



(2) 



where <j>(r, A) is the SSS solution of the one-component 
scalar field associated with the "mean-square" scalar 
charge A = ^/^ a (A a ) 2 , and the x' a s are integration 
constants. The energy associated with these solutions, 
when finite, is degenerate on spheres of radius A in the 
N-dimensional A a space and equals the energy of the as- 
sociated one-component scalar field with scalar charge A. 
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Moreover, the admissible scalar field models with soliton 
solutions have been classified into six families with regard 
to the behaviors of the SSS field strengths at the center 
(r = 0) and as r — > oo. Such families arc the combina- 
tions of five characteristic behaviors. Classes A-l and A-2 
correspond to infinite (but integrable) and finite soliton 
field strengths at the center, respectively. Classes B-l, 
B-2 and B-3 correspond to asymptotic dampings of the 
soliton field strengths which can be slower than coulom- 
bian, coulombian or faster than coulombian, respectively. 
For every class we obtained the associated behavior of 
the lagrangian function f(X) around the corresponding 
value of X. Finally, all the scalar soliton solutions of ad- 
missible models were shown to be linearly stable against 
vanishing scalar charge perturbations. 

In the present paper we shall determine the class 
of generalized gauge-invariant models supporting sim- 
ilar electrostatic spherically symmetric (ESS) soliton 
solutions. Let us begin with the simpler case of 
electromagnetic fields and go next to the analysis of non- 
abelian gauge fields. In the former case we write the 
general expression for a relativistic gauge-invariant la- 
grangian density as 



L = <p{X,Y), 



(3) 



where ip{X, Y) is a given function of the two field invari- 
ants 



X = -If^F"" = E 2 - H 2 
Y = -If^F*^ = 2$. it, 



(4) 



which is assumed to be continuous and derivable in the 
domain of definition (f2) of the X — Y plane (5ft 2 ). We 
also assume the condition ip(X,Y) = ip(X,—Y) to be 
satisfied, in order to preserve parity invariance. As in 
reference [1] we shall call "class- 1" field theories the mod- 
els for which ip(X, Y) is defined and regular everywhere 
(Q = 5ft 2 ) and "class-2" field theories those with f2 C 5ft 2 , 
(0, 0) G 51 and fl open and connected. Other models are 
excluded from this analysis. The original BI electrody- 
namics is a class-2 theory corresponding to the particular 
choice 



Vbi{X,Y) 



1 - n 2 X - ^-Y 2 



(5) 



which reduces to the Maxwell lagrangian <p(X, Y) = 
in the low-energy limit (/i — ► 0). We look now for the 
conditions to be imposed on the functions tp(X, Y) in 
order to define consistent general electromagnetic field 
theories. First of all let us consider the energy density 
obtained from the symmetric (gauge-invariant) energy- 



momentum tensor: 



■ dLP :E 2 + 2^E.H-<p(X,Y) 



dX' 



, 9<P 
J dY J 



The requirements of the positive definitencss of the en- 
ergy and the vanishing of the vacuum energy lead to the 
following necessary conditions to be imposed on tp(X, Y) 



tp(0, 0) = ; <p{X, 0) < (VpT < 0, Y = 0) e Q) 
dip 



^>o (V(x,Y)efi), 



(7) 



aside from the minimal necessary and sufficient condition 
(to be satisfied for V(A, Y) e ft) 

p s > (X + Jx^TT 2 )^ +Y^- <p(X, Y) > 0. (8) 

For the purposes of this study we shall consider "admis- 
sible" only the models satisfying these conditions. Then 
the admissible models are the solutions of the linear in- 
homogeneous partial differential equations obtained by 
equating the expression in Eq.® to any positive func- 
tion and satisfying the conditions ([7]). 

The field equations obtained from Eq. ([3]) take the form 



dip 



°^ OX 



dip 
dY J 



0. 



(9) 



For central electrostatic fields, Y = 0, and these equa- 
tions have a first integral, which may be written in 
terms of the electrostatic potential Ao(r) (A = 0, E(r) — 
-VA (r)) as 



r2 dAo£^ {XY 



dr dX 



0)=?, 



(10) 



where X — (^r) 2 = E 2 {r) and q is the integration con- 
stant. The solutions of this equation, when substituted 
in Eq.Q, lead to a Dirac S distribution of weight Anq. 
The definition of the electric charge associated to a given 
field 



1 f ( dip - dip 

^j d ^\dX E+ dY^ 



H 



(11) 



allows to identify this constant as the electric charge of 
the solution [l[. With the identification cf>(r) = A (r), 
Eq. pT)]) coincides with the SSS field equation of a one- 
component scalar field model (Eq.([T]) without indices) 
with a lagrangian density given by 



L scalar = f(X) = -ip(-X, Y = 0). 



(12) 



Conversely, we may associate to every admissible scalar 
model defined by a lagrangian density f(X), a family 
of admissible electromagnetic field models defined by la- 
grangian densities ip(X, Y) which satisfy Eqs. (fT2)l and the 
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conditions |7J) and {8]| . The absolute values of the elec- 
trostatic central field solutions of all these generalizations 
(\E(r,q)\) have the same form, as functions of r, as the 
SSS field solutions (<f> (r, A)) of the original scalar model. 
In calculating the electrostatic energy, by integration of 
Eq.©, and comparing to the energy of the associated 
scalar solitons with the same integration constant we are 
lead to the relations 

e e {q) = 8wq [A (oo, q) - A o (0, q)} - e s (q) = 2e s (q), (13) 

(indices e and s stand for electric and scalar fields, re- 
spectively). We conclude that the analysis of the con- 
ditions determining the scalar field models with finite- 
energy SSS solutions, performed in Ref.[8], allows the 
determination of families of admissible electromagnetic 
field models supporting finite-energy ESS solutions, and 
the energy of such electrostatic solutions of charge q is 
twofold the energy of the corresponding scalar soliton 
with scalar charge A — q. Equivalently, the scaling of 
energies in the scalar case (e s (A) = A 3 / 2 e s (A = 1)) leads 
to the relation q = (2) 2 / 3 A between the charges of an 
electrostatic soliton (q) and the associated scalar soli- 
ton (A) of equal-energy. Moreover, the classification of 
the scalar models through the central and asymptotic 
behaviors of their soliton solutions, induces a similar 
classification of the gauge models through the behavior 
of their finite-energy ESS solutions. The conditions on 
the lagrangian functions ip(X, Y) leading to the differ- 
ent classes of finite-energy solutions are immediately ob- 
tained, through Eq. (fT2|) . from those leading to the clas- 
sification in the scalar case. 

However, although in the scalar case the finite-energy 
SSS solutions of admissible models are always linearly 
stable, this is not so for the generalized gauge-invariant 
models. These models must satisfy supplementary con- 
ditions in order to support ESS soliton solutions. To find 
these conditions let us analyze the linear stability of the 
ESS solutions, which requires their energy to be a min- 
imum against charge-preserving perturbations. Let us 
consider a finite-energy ESS solution of the field equa- 
tions © (Eq (r) , Hq — 0) and introduce a small regular 
perturbing field (Ei(f,t), Hi(r,t)) with vanishing elec- 
tric charge density. To the first order, the perturbing 
fields satisfy the equation 



V.S = 0, 



where 



£ = 



dip 



Ei 



'dX* 



(Eo-Ei)Eq, 



(14) 



(15) 



and the index in the derivatives means that they are 
calculated for the unperturbed solution. The first-order 
variation of the energy functional under these small per- 
turbations, obtained from the integral of Eq.©, becomes 

Are = -2 / d 3 rV. \a°£] + 2 [ d 3 rA°f.T,, (16) 



where we have introduced the time- like component of the 
vector potential for the solution (Eq = — S7A°). From 
the linearized field equations (fT4|) and the asymptotic 
behavior of £ (which results from the integration over all 
space of the same Eq. (fT4"]) ) , we see that the first variation 
of the energy functional vanishes, which is an extremum 
condition. The expression of the second-order variation 
of the energy functional, obtained from © together with 
the expansion of the field equations up to the second 
order, reads 



A 2 e = 



+ 



tar 



dX 



+ 



(17) 



By deriving Eq. (|10[) with respect to r and taking into ac- 
count the monotonic character of |i?o( r )l we can show the 
positivity of the first integral for any admissible model. 
Thus the requirement 



dX dYf 



(18) 



to be satisfied in the range of values of X (Y = 0) de- 
fined by the ESS solutions, is a necessary and sufficient 
condition for the linear stability of the finite-energy ESS 
solutions of the admissible models. This stability cri- 
terion goes beyond the widely used Derrick's necessary 
conditions [111 ]. 

We shall now proceed to extend the analysis to 
non-abelian gauge field theories of compact semisimple 
Lie groups of dimension N. We use the following con- 
ventions for the tensor field strength components in the 
algebra and their duals, in terms of the gauge fields A ayL 
and the structure constants C a bc- 



Eajiv — dp,A av — d u A a ^ g ^ ' CabcA^^Ac 
1 



be 



(19) 



whose components define the fields E a , H a in the usual 
way. In building the lagrangian density functions in this 
case there is an ambiguity in the calculation of the traces 
over the matrix- valued fields [l(|. We adopt here the 
ordinary definition of the trace and introduce the field 
invariants X and Y as 



a a 
y=-\ J2( F ^Kn - 2^ (Ea-Ha) , (20) 



where 1 < a < N. The lagrangian density is now defined 
as a given continuous and derivable function <p(X, Y) 
which, for parity invariance, is assumed to be symmet- 
ric in the second argument (tp(X, Y) — ip(X,—Y)). We 
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also introduce the distinction between class- 1 and class-2 
models with the same criteria as in the electromagnetic 
case. The associated field equations read now 



E 

b 



D, 



abfi 



9( P T?^ 

ox » 



dY b 



0. 



(21) 



where D abll = 5 ab d^, + gJ2 c C ab c A Cfl is the gauge covari- 
ant derivative. The energy density associated to the sym- 
metric energy-momentum tensor is 



T n s , 



(22) 

As can be easily verified, the conditions to be imposed on 
the lagrangian density for "admissible" models (positive 
definitcness of the energy density and vanishing vacuum 
energy) take the same form (0, ([8])) as in the abelian 
case. 

Let us consider the electrostatic spherically symmetric 
solutions of the field equations (f2~Tj) . We look for fields of 
the form 



K(r) = -V(0 Q (r)) = -0 o (r)- 



0, (23) 



where the functions </> a (r) = A®(r) can be identified as 
the time-like components of the gauge potential in the 
Lorentz gauge. When substituted in the field equations 
(12T1) we are lead to 







be 



CabcMrWAr) = 0. 



(24) 



The first group of equations leads to the set of first inte- 
grals 



2 dl P 



(25) 



where X = ^2 c <fic( r ) and Qa are the integration con- 
stants, which must be interpreted as the color charges 
associated to the soliton. Indeed, as in the abelian case 
the color charges associated to a given field read 

which now include the external source charges and the 
charges carried by the field itself. The latter ones vanish 
for the ESS solutions and it remains only the former, 
associated to Dirac distributions of weight 47rQ a resulting 
from the substitution of Eqs.(j25|) in the first of (|24|) . 

Equations (|25|) coincide with the field equations ([1]) 
for a multicomponent SSS scalar field, whose lagrangian 
density is given by 

L = f (j^d^d^ = f(X) = -<p(-X,Y = 0). 

(27) 



Consequently their solutions coincide with those of 
Eq.([T]) for the associated N-components scalar field the- 
ory, and reduce finally to the solution of the associated 
one-component scalar problem. These solutions are then 
given by Eq.® with the substitution A a — » Q a . They 
must also satisfy the second set of equations (|24p . This 
restricts the possible values of the constants Xa and leads 
to the following final expression of the general ESS solu- 
tions 



0a(r) = ^f(0(r,Q)+ X ) 



(28) 



where Q — ^Y^ a (Qa) * s tne "mean-square" color charge 
and x is an additive constant which is the same for all 
the components of the potential. In terms of the fields 
these solutions read 

E a (r) = \r,Q)~ ; H a = 0. (29) 

Q r 

The associated energy is obtained from the integration 
of Eq . (f22|) . It is now straightforward to show that (as in 
the many-components scalar case) this energy, if finite, 
is degenerate on spheres of radius Q in the color charge 
space and is related to the energy of the corresponding- 
scalar soliton of charge A = Q through the same equa- 
tion (fl"3]) of the abelian case. Moreover, the criteria ob- 
tained in Rcf.[8| determining the classes of admissible 
scalar field models with SSS soliton solutions determine 
also, through Eq. ((27|) . the classes of admissible general- 
ized non-abelian gauge models supporting finite-energy 
electrostatic particle-like solutions. 

As in the electromagnetic case we must look now for 
the conditions leading to the stability of these ESS solu- 
tions. Nevertheless, as already mentioned, a new qualita- 
tive effect arises in the present situation, which is related 
to the intrinsic color charge density associated to the 
fields perturbing the ESS solutions. By defining E la (r, t), 
Hi a (r,t) and A^ a (r,t) as the perturbing fields and po- 
tentials respectively, we must require that the total color 
charges of the perturbed solution remain unchanged. To 
lowest order the first of Eqs. ([M)l leads to 



V.E Q = Ain Q , 



(30) 



where the corrections to the color charge density distri- 
butions introduced by the perturbations, Ain a (r, t), and 
the term E a are given by 

dip 



Ain a 



dX 



C a bcAib-E c 



be 



<>f Eia + 2^^{E b .E lb )E a . (31) 



dX 



These corrections must satisfy the charge-preserving con- 
dition j- J d^rAiria — 0, which determines the asymp- 
totic behavior of E a through the integration of Eq. (|30|) . 
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The lowest-order perturbation of the energy, obtained 
from the integral of ([22")) . using (|30|) and this charge- 
preserving condition becomes 



Aie 



(32) 



which must be interpreted as the interaction energy of 
the field of the solution with the first-order correction 
to the charge density distribution introduced by the per- 
turbations. This contribution vanishes identically due to 
the second set of equations ([24]) . Consequently, the first 
variation (|32j) vanishes, and the solution is an extremum 
of the energy functional. 

In the same way we obtain the second variation of the 
energy by perturbing up to the second order the energy 
functional and the field equations. The final expression 



A 2 £ 



(far 



(far 



Yl ^la + 2 ( ^Oa -Ela 



dX, 



dip 



X] Hi* ~ 2 ( X Eoa-Hia j ( 33) 

n \ n I 



The first integral is always positive. Thus the require- 
ment of positivity of A 2 £ leads to the condition 



8X > dY*- 



(34) 



to be satisfied in the domain of values of X (Y = 0) de- 
fined by the solutions. This condition, which is formally 
the same as in the abelian case, is necessary and suffi- 
cient for the linear stability of finite-energy ESS solutions 
of admissible non- abelian gauge models. 

As already discussed in Ref. the linear stability of 
these non-topological solitons does not guarantee their 
general stability, intended as the preservation of the "soli- 
ton identity" when interacting with strong "external" 
fields. Only in some few special models with topologi- 
cal soliton solutions (mainly in one-space dimension) the 
existence of conserved topological charges allows the "de- 
tection" of the presence of solitons in a field configura- 
tion [l^ . But despite some tentative approaches to this 
question [13j, the development of general methods for 
the analysis of the stability and dynamical evolution of 
strongly-interacting solitons in three space dimensions re- 
mains an unsolved problem. Nevertheless, for low-energy 
interactions (weak external fields), the identity of the 
soliton is preserved if (and only if) linear stability holds. 
In this case an expansion in the characteristic energy 
scale of interaction (or, equivalently, in the character- 
istic value of the intensity of the external field) becomes 
possible. At the leading order this expansion allows the 
interpretation of the field dynamics in terms of particle- 
field force laws and radiative behavior , ■ 



Let us conclude with a brief comment about the impor- 
tant question of the radiative solutions of the generalized 
gauge field theories. In this work we have focused on 
the static soliton solutions of these theories and, from 
this point of view, the Born-Infcld-like models are only 
particular examples of a large family, without special dis- 
tinguishing properties. However, the BI original abelian 
version has very special properties from the point of view 
of wave propagation. Indeed, as can be easily verified, all 
these theories exhibit plane wave solutions propagating 
with the speed of light but, owing to the non-linearity, 
there are also other radiative solutions propagating with 
more complex dispersion relations. In most cases the 
energy of such waves tends to cumulate, evolving to- 
wards spatially-singular configurations, with formation of 
shocks after some critical time. A system of field equa- 
tions whose wave solutions are free from this evolution is 
called "completely exceptional" . In analyzing this ques- 
tion for the case ofgeneralized electromagnetic field the- 
ories, G. Boillat 14| established that, from among the 
family of all admissible models belonging to the class B- 
2 defined here (asymptotic coulombian behavior of the 
ESS solutions; the condition §$(X = 0, Y = 0) = 1 im- 
posed on the lagrangian excludes the classes B-l and B-3 
from the Boillat analysis) only the Born-Infeld model is 
completely exceptional. However, the analysis of mod- 
ern versions of the Born-Infeld model, in the contexts of 
Kaluza-Klein or string theories, concludes that this ex- 
ceptionality is lost in the generalizations of the classical 
version [la ]. 

In summary, by extending the methods developed in 
Ref. [|| for the scalar field models, the present analysis 
allows the explicit determination of any admissible 
generalized gauge field model supporting ESS soliton 
solutions. More details, examples and applications of 
these results will be published soon [lij ]. 
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